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In this paper we introduce a new method to design control laws for non-linear un- 
deractuated systems. Our method will often work in situations where standard control 
techniques fail. Even when standard techniques apply, we believe that our approach will 
achieve better performance. This is because we produce an infinite dimensional family of 
control laws, whereas most control techniques only produce a finite dimensional family of 
control laws. We will describe the problem and our solution invariantly, using differential 
geometry and in local coordinates. We hope that this paper will be useful for both math- 
ematicians and engineers. We include an abstract example of a system which is open-loop 
unstable and cannot be stabilized using any linear control law, and demonstrate that our 
method produces a stabilizing control law. We also apply our method to the inverted 
pendulum cart to compare it with some standard control techniques. 

An important problem in control theory is how to modify a system of ordinary dif- 
ferential equations so that solutions to the equations satisfy some property. In order to 
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describe the class of differential equations that we consider, let Q denote the configu- 
ration space. The configuration space is a finite dimensional manifold that represents, 
for example, every possible position of a mechanical system. Let g G T(T*Q ® T*Q) 
be a metric. This is the mass matrix in a mechanical system. Let c, / : TQ — > TQ be 
fiber-preserving maps. These maps are not assumed to be linear, but we do assume that 
c is odd, i.e., c(— X) = — c(X). In a mechanical system c will represent the inherent 
dissipation and / will represent the input forces. Finally, let V : Q — > IL This will 
represent the potential energy of a mechanical system. The differential equation that we 
consider in this paper is 

V^7 + c(t)+ grad 1 V = f(j). (1) 

Let P G F(T*Q <S> TQ) be a ^-orthogonal projection. We consider the situation where a 
constaint P(f) = is imposed. A system is called underactuated if P 7^ 0. A mechanical 
system is underactuated if we require that control forces are zero in certain directions. 
Many problems in chemical, electrical and mechanical engineering may be formulated 
as follows: find a function / with P(f) = so that solutions to Equation (1) have 
some specified property. For example, solutions with initial conditions in some region 
will remain close to some path, or as a different example, some point in TQ will be an 
asymptotically stable equilibrium. 

In local coordinates x = (x , . . . , x n ) Equation (1) reads 

x k + T k j (x)x i x j + c k (x,x) +g tk ^~ = f k (x,x), (2) 

where Y^Ax) are the Christoffel symbols (of the Levi-Civita connection) associated to 
the metric g, [Hicks, 1965], and g lk is the inverse matrix to gj-j- In Equation (2) and 
throughout the rest of this paper we are using the summation convention, that repeated 
indices are summed from 1 to n. As we have mentioned, we assume that c k (x, — x) = 
—c k (x, x) and the projection of f k to some specified set of directions must be zero. 

The main question addressed in this paper is how to find a function, /, with P(f) = 0, 
so that solutions to Equation (1) satisfy some preassigned conditions. Our approach to 
this question is to find functions g, c, V and / so that solutions to Equation (1) are 
automatically solutions to 

V^7 + 0(7) + grad^V = 0. 
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This will clearly be the case if 

/(X) = V X X - V X X + grad^V - gT^d^V + c{X) - c{X), (3) 

for every vector field X. The condition P(f) = then becomes a system of nonlinear 
partial differential equations for g, c, and V. Notice that constant multiples of g, c, and 

V satisfy P(f) = even when P has full rank. Thus, one would expect many solutions 
when P does not have full rank. Separating P(f) = into terms which are quadratic in 
the velocity, independent of the velocity or odd functions of the velocity gives 

P(V x X-V x X) = 0, (4.1) 

P(grad 7 V - grad^V) = 0, (4.2) 

P(c(X)-c(X)) = 0. (4.3) 

We will look for solutions to these matching equations with g non-degenerate so that 
g(X, Y) = g{\X, Y) with A G T(T*Q <g> TQ). It is clear that A has to be g self-adjoint, 
i.e., g(XX, Y) = g(X, XY). We will derive a linear system of partial differential equations 
for A which must be satisfied if g is to solve Equation (4.1). 

To derive the equations for A we will use the relation between the connection, Lie 
bracket (commutator) and metric. In fact, we only derive equations for A| Im p . 

It is known, [Hicks], that the covariant derivative V compatible with the metric g is 
determined uniquely as a bilinear operator which associates to any pair of vectors X and 

Y a third vector V xY so that the following equations are satisfied: 

X(g(Y, Z)) = g(V x Y, Z) + g(Y, V X Z) , (5) 

V x Y-VyX=[X,Y], (6) 

where [X, Y] is the Lie bracket (commutator) of X and Y. Using the above properties 
of the covariant derivative we get: 

2g( V X Y, Z) =Xg(Y, Z) + Yg(Z, X) - Zg(X, Y) 

+ g([X,Y],Z) + g([Z, X] , Y) - g([Y, Z],X). (7) 



This may be solved for Vx^ since g is non-degenerate. 
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Proposition 1. If g = gX and g satisfies the matching condition P(Vx^ — Vx^) = 0, 
then X satisfies: 

V#A| Im pm = 0, (8) 

or, equivalently, 

g(VzXPX,PX)-g(\PX,V z PX)=0. (9) 
Proof. We begin by polarizing the matching equation 

P(V X Y - W X Y) = l -P{V x Y + Vyl - V X Y - V Y X) 



1 
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P(V x +y(X + Y)- V x+y (X + Y)) 
-P(V X X - V X X) - P(V Y Y - V Y Y) 



0. 



The first line is true because the covariant derivatives are torsion free (Equation (6)). The 
second line is true because the covariant derivatives are bilinear (Equation (7)). Now, 

= 2g(P(V XPX Z - V X pxZ),X) 
= 2g(V XP xZ - VxpxZ, PX) 
= 2g(V X pxZ, PX) - 2g(V X pxZ, XPX) 
= XPXg(Z, PX) + Zg(PX, XPX) - PXg(XPX, Z) 

+ g([XPX, Z],PX) + g([PX, XPX],Z) - g([Z, PX],XPX) 

- XPXg(Z, XPX) - Zg(XPX, XPX) + XPXg(XPX, Z) 

- g([XPX, Z],XPX) - g([XPX, XPX],Z) + g([Z, XPX], XPX) 
= XPXg(Z, PX) - PXg(XPX, Z) + g{[PX, XPX],Z) 

- g([Z, P X], XPX) + g([Z, XPX],PX) 
= g(VxpxZ, PX) + g(Z, VxpxPX) 

- g(XPX, V PX Z) - g(Z, V PX XPX) 

+ g([Z, XPX],PX) - g([PX, XPX],Z) - g([Z, PX],XPX) 
= g(V z (XPX),PX) - g(V z PX, XPX) 
= g((V z X)(PX),(PX)) 
= (V z gX)(PX,PX). 



The covariant derivative is extended to all tensors by requiring every reasonable prod- 
uct rule to hold. For example, Z(X(X)) = V Z (X(X)) = (V Z X)(X) + A(V ' Z X) and 
X(g(Y,Z)) = W x (g(Y,Z)) = (V x g)(Y, Z) + g(V X Y, Z) + g(Y,V X Z). We used these 
two product rules in the last two lines and Equations (5), (6) and (7) in the previous 
lines. D 

Now that we have equations for A| Im p , we will derive equations for g in terms of 

A llm P- 

Proposition 2. Ifg = 'gX, then 

XPXg(Z, Z) + 2g([Z, XPX],Z) = 2Zg(PX, Z) - 2g(PX, V Z Z) (10) 

Proof. 

XPXg(Z,Z) + 2g([Z,XPX],Z) 

= 2 g(V XPX Z, Z) + 2g( [Z, XPX] , Z) 
= 2g(V z XPX, Z) 

= 2Zg(PX,Z) -2g{PX, V Z Z). □ 

The matching equation for the potential energy may be expressed as a linear partial 
differential equation using A. 

Proposition 3. If the matching Equations (4) are satisfied and g = gX, then 

XPXiV) = PX(V). (11) 

Proof. XPXiV) = (dV)(XPX) = gQrad^V, XPX) = g^grad^V^X) 

= g(P gToTd^V.X) = g(P grad^X) = PX(V). □ 

By solving the matching equations, we find every control law that will result in dy- 
namical behavior equivalent to a system of the form: 

V^7 + c(7) + grlTd^V = 0. 



This is a very large collection of vector fields on TQ. Thus far, we have shown that 

every solution to the matching equations may be found by first solving one system of 
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linear partial differential equations and then solving a different system of linear partial 
differential equations. 

The previous three propositions suggest a natural approach for solving the matching 
equations. First, solve Equation (9) for XPX, then solve Equation (10) for g, and Equa- 
tion (11) for V. Finally solve the algebraic Equation (4.3) for c. Every solution to the 
matching equations can be found in this way. What is not clear is whether every set of 
functions generated in this way is a solution to the matching equations. The problem is 
that once we have a solution to the A-equation, (Equation (9)), and then a solution to 
Equation (10), it is not clear that A can be extended from Im P so that the condition 
g = gX holds. We will address this question for general P in a future paper. We will next 
show that g = 'gX holds if it is true at the initial conditions, provided that the rank of P 
is one. 

Proposition 4. Let X be a non-zero vector field in an open set U, generating Im P , and 
let XPX be a non-vanishing solution of the X Equation (9) in U. Let E be a codimension 
one non- characteristic hypersurface in U . Assume that g is defined, non- degenerate, and 
satisfies TjXPX = gPX on E, and V is defined on E. Then there is a unique solution to 
the matching Equations (4-1) o,nd (4-2), and a unique extension of X away from Im P in 
a neighborhood o/E. Furthermore, gX = g everywhere in the neighborhood. 

Proof. Notice that E is non-characteristic for the ^-Equation (10) if it is non-characteristic 
for the l^-Equation (11), and visa versa. The functions g and V may be found in a neigh- 
borhood of E using the method of characteristics. The main content of this proposition is 
that gPX = 'gXPX everywhere in the neighborhood. After polarization, the (/-Equation 
(10) reads: 

XPXg(XPX, Z) + g(XPX, [Z, XPX}) 

= Zg(XPX, PX) + XPXg(PX, Z) - g(V z XPX, PX) - g{V XPX Z, PX) 
= XPX g(PX, Z) + g{PX, [Z, XPX}) + g(XPX, V Z PX) - g(V z XPX, PX). 

This last line follows after applying the product rule to the first term of the previous 
line and the torsion-free condition (6) to the last term of the previous line. In view of 
the A-Equation (9), it is clear that g^XPX = gPX is a solution of this equation. This 
solution is unique in a neighborhood of E. Finally, since ^ is non-degenerate on E, it is 
non-degenerate in a neighborhood of E, so we may extend A to the full tangent space by 
A=(?)-V □ 

6 



EXAMPLE 

In order to demonstrate and clarify our method, we will apply it to an inverted pen- 
dulum cart as an example. 




y 

Figure 1 
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The (mass) metric for the cart system depicted in Figure 1 is: 



g = (M + m)dy 2 + 2m£ cos(9)d9dy + (m£ 2 + I)d6 2 , 
and the potential energy is: 

V = mg£cos(9). 

Here 

M is the mass of the base of the cart, 
m is the mass of the pendulum, 

£ is the length from hinge to the center of mass of the pendulum, 
/ is the moment of inertial about the center of mass, and 
g is the acceleration due to gravity. 

By a change of length, time and mass scales, the metric and the potential energy may 
be transformed into 

g = dx 2 + 2bcos(9)d9dx + d9 2 , V = cos(9), (12) 

where b = m£{M + m)~^{m£ 2 + E (0, 1) is a dimensionless parameter. 

The problem is to find the force applied to the base of the cart as a function of the state 
of the cart so that the origin (9 = 0, x = 0) will be an asymptotically stable equilibrium 
of the resulting dynamical system. 

In this application, we may only directly apply force in the x direction, so P is pro- 
jection onto the ^ direction, i.e., 

d 

P= (bcos(6) dx + dO)®—. 

ov 

Projecting system (1) onto the coordinate directions gives: 

P(V^7 + c(7)+ grad 1 V) = 0, 

(0 \ ( d \ (13) 

— ,V^ 7 + c( 7 )+ grad^Vj =u = g( — ,f) . 

We will assume that there is no inherent dissapation: c = 0. Then, in local coordinates, 
the equations read: 

9 + bcos(6)x - sin(0) = 

•• (14) 
bcos(9)9 + x - bsm(9)9 2 = u. 



Equation (3) expresses our solution to this problem in terms of functions, gij, V and 
Cj, i,j = 1 or 2, which satisfy the Equation (4). In local coordinates, these matching 
equations read: 



(#22 - 6cos(6>)#i 2 ) 
(#22 - bcos(6)g 12 ) 



8g 



li 



d6> 

dx 



+ (6cos(0)<7u -#12) 



,<9#12 <9#11 



d6> 



5a; 



/22 



<9# 



+ (&cos((9)^n -#12)- 



/22 



<9x 



= 0, 

= 0, 



(#22 - 6cos(6')^i2) :: ^ + (6cos(6>)#n 



^12)^ = 



<9U 8V 
(#22 - bcos(d)g 12 ) — + (6cos(0)(7ii -#12)-^ = - sin(0)(<7u<7 2 2 - (#i2) 2 ), 



and 



ci + 6cos(6>)c2 = 0. 

We have slightly simplified these equations by multiplying by the determinant of g. (The 
terms on the right hand side of each equation would be non-zero if the cart was on a 
hill, and the hinge was rusty. Therefore, multiplying by the determinant of g would not 
be a great simplification.) These equations are messy, but the same functions may be 
found by solving the equations in Propositions 1, 2, and 3 without ever writing out the 
matching equations. Take X = in Proposition 1 (so PX = J^) and let 
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89 



d 



a — = a— + flor- 



ae 



dx 



Using the bilinearity of the covariant derivative and the product rule, the A-Equation (9) 
becomes: 
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d d \ da 



and 



<96>' dd 89 



d d \ da 



+ # 



d d \ d\i 



00' 89 dx 



dx 1 d9 d9 



d d \ dn 



+ 9 V a 



d d 



dx d9 J dx 



™ dx' d9 



d d 



V- 9 



— V — 

dx' ^89 



H = 



<7h^>^h^ + <M^,™hf + <MVa— ,— )f l -g[— ,Va — U = 



^8x : 89 
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dx' ^d9 



Using Equation (7), it follows that: 

da 



d\i 



QQ + &cos(0)-^ + &sin(0)/z = O, 

da . /n ,da 
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By the Frobenius Theorem, or equivalently checking that mixed derivatives are equal, we 
see that this system has a solution if and only if b sin(0) |^ = 0. The second equation then 
implies that a only depends on 9, so the first equation reduces to an ordinary differential 
equation. The general solution to this system of equations is, therefore, 

a = a(9), n = /io cos(0) — - cos(#) / — sec 2 6 d6. 

b J o9 

For the rest of this example, we will use the particular solution: 

a = (To, jtt = /xo cos(#) . (15) 

It is easiest to find g using the basis J|, A (-§$)■ By the relation #A = g, we already 
know g(Hm)>m) and ? (A (^) A (^)), so we only need to find g(-§g,-§g). The Lie 
bracket, 



^ aA 
d9 J K d9 ) 



d d d 

= -/x o sm(0) — = a tan(6')— - tan(0)A I ^ 



may be plugged into the ^-Equation (10) to get: 

a o ^ + A»o cos(9) + 2a tan(0)£n - 2 tan(0) = 0. 
o9 ox 

The surface 6* = is non-characteristic, so the above equation may be solved with initial 
data gn(0,x) = h(x). 

The flow of the vector field XPX is given by: 

6 ' = (To, x = /xo cos(0), 

or 

9 = dot, x = a H sin(<T i), 

to 

using 6> = at t = 0. With this flow, the # equation becomes: 

(it 

so 



+ 2<To tan(<Tot)^n — 2 tan(<ro£) = 0, 



gu = —+ (h(a) - —) cos 2 (a t) = — + ( h ( x - ^ sm(9)] - — ) cos 2 (#). 
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Again, we only take a particular solution: 

9u = — + rcos 2 (6»), 
where r is a constant. The equations 

bcos(9) =g{^^)=9 (A(^), ^) = *09i2 + fJ'O cos(0)& 2 

and 

1 = 9 (w We) = g ( A( l } ' We) = aodl1 + ^ cos{e)dl2 

may now be solved for the remaining terms of g: 

912 = rcos(6>), g 2 2 = h -g?". 

We can compute in the same way. Using the XX P flow and initial data V(0, x) 
w(x), the F-Equation (11) may be written as an ordinary differential equation: 

dV . . , 
— = -smM). 



Thus, 



i> = J_( CO s((jot) - 1) + w(a) = — (cos(0) - 1) + w ( x - — sin# 
O"0 o"o V °"o 



Finally, solving for c, we get: 

bcos(e)—- — 

with an arbitrary function K. 

The functions 'g, V and c~ produce a family of control laws via equation (3). In this 

example, the control laws depend on four functions, ctq(^), h(x), w(x), and K(0,x,e,x), 

and one constant, hq. We have already chosen ao and /i to be constants. We will address 

the question of how to choose the unknown functions and parameters in order to best 

meet specific design criteria in a future paper. For now, we will just pick elementary 

functions that will insure that 6* = 0, x = is an asymptotically stable equilibrium. If 

ll 



(0,0) is an equilibrium, gradV(0,0) must be zero, so w'(0) must be zero. One standard 
way to insure that (0, 0) will be an asymptotically stable equilibrium is to pick positive 
definite g : V and gc, then H(X) = ^g(X, X) + V will be a Lyapunov function with time 
rate of change — g(c(X),X). The Hessian of V at (0,0) is 



n 2 - 

^(o,o) - 



4^(0) - - 



w (0) 



So we should require: 



- ^V'(o)^ 

w"(0) 
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w"(0) > and det Df 0fi) V = -—w"(0) > 0. 



We will choose cr < 0, and w(z) = \w\z 2 with w\ > 0. 

The mass matrix should also be positive definite. Thus we should have gn > and 
#22 > so that r > and + o"o r > 0. In addition, the determinant of the mass matrix 
will also be positive: 

^ + ^cos 2 W + ^r>0. 
Rearranging and using the fact that a < gives: 

— aobfior cos 2 (9) > a 2 r + b^Q. 

So, 



0-5 r + bfi 



and 



cos 2 (0) > 



aobr cos 2 (9) 
a 2 r + bjxo 



>0, 



(16) 



-a bfx r 

Assuming /io > 0, r > and 00 < 0, the mass matrix will be positive definite provided 
condition (16) holds. When the mass matrix and potential energy of the model system 
are positive definite, the point (0, 0) will be Lyapunov stable for any positive semi-definite 
go. In particular, when K = (i.e., c = 0), we will get Lyapunov stability. We call 
controllers with c = conservative controllers. A direct computation shows that: 



g(c(x),x) 



Kg 
K 



/m 9 d-d .8 



(bcos(9)g 11 - gi 2 )9 + (bcos(9)g 12 - #22)^] 
a x). 



= (det g) ■ K ■ (fi cos 99 
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This will be positive semi-definite if K = $(#, x)(jj, cos 99 — a x), where $ is any pos- 
itive function. (We can do no better because XPX = a -^ + hqcos9-^ will be a zero 
mode of any admissible c.) We will see that this is sufficient to prove that (0, 0) is an 
asymptotically stable equilibrium of the controlled system. Equation (3), equation (7) 
and our values for g, V and c combine together to give: 



"'' <K ' ' j ( x- —sm(9)] + detg$(9,x) {n Q cos(9)9 - a x), 



a detg \ a 



(17) 



where 



det# = l-& 2 cos 2 (6>), detg = — + — cos 2 (#) + (18) 

(T /i IH) 

Thus, we have almost proved the following result. 

Proposition 5. Let constants /j,q, (Jq, r , and wi satisfy the conditions 

a 2 r + bun 

^o>0, (7 <0, Wl >0, 1 > , W . 19 

-a b^i Q r 

Let <&(9,x) be any strictly positive function. 

Then (0, 0) is an asymptotically stable equilibrium of the controlled system (14) with 
the control law defined by (17), (18). 

Proof. Define the controlled Hamiltonian H = ^g(j, 7) + V. We have shown above that 
the Hessian of H is positive definite in some neighborhood of (0, 0) when conditions 
(19) hold. The time derivative of H is: 



dV ' dt \2' 

= g(^ii+ gradV,^) 
= -<7(c(t),7) 

= -detg$>(9,x) (fi cos(9)9 - a x) 2 . 

Check that there is no solution to the controlled equations satisfying /j,q cos(9)9 — o-qx = 0. 
Since $ > by assumption and detg|^^ > for all (9,x) sufficiently close to (0,0), 
the Hamiltonian H decreases along the solutions of (17), and, therefore, may serve as 
a Lypunov function. The (local) asymptotic stability follows from Lyapunov's Theorem. 
□ 
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Remark 1. It is not possible to construct a control law so that (0, 0) is a globally asymp- 
totically stable equilibrium. 

Indeed, the solutions of the controlled system for such a control law would produce 
a continuous function, F : [0, oo] x TQ — > TQ so that Fo is the identity map and 
is the constant map. In other words, the flow of a vector field with a globally 
asymptotically stable equilibrium is a contraction. However, for the inverted pendulum 
cart, TQ ~ 5" 1 x IR 3 , which is not contractible. 

Since the inverted pendulum cart cannot be globally stabilized by a control law, we 
can only try to maximize the size of the basin of attraction. We will compare a special 
case of our nonlinear control law with a linear control law that has been implemented on 
an inverted pendulum cart in our lab. Our cart has 

M = 5.02 Kg m= .454 Kg 
t = .425 m 1 = .11 Kg m 2 

giving b = .188. The linear control law that we obtained by placing poles at —5, —6, and 
a double pole at —2 is: 

uiin = 10210 + 115.8x + 918.50 + 158. 2x . 

This is only an approximation of the actual control law that is running in the lab. In the 
lab we had to model the DC motor, the observer, and use a discrete control law. 

Thus far, we have chosen each of the arbitrary functions in our non-linear control 
law to be a constant, and we took w(x) = \vj\x 1 , since it has to have positive second 
derivative. As a rough guess we took fj, = 10, oo = —.05, r = 1000, and w\ = 1.5. 
We also set the function $ = 1. These constants satisfy conditions (19). We chose the 
functions to be constants when possible just to simplify the exposition. 

Remark 2. This would not be the best choice for engineering applications. In particular, 

— oo controls the coefficient of x in our control law. In order to stabilize the system after 

a large angular disturbance, we would like —o~q to be small so that the cart will be free 

to accelerate under the pendulum. On the other hand, —ao needs to be much bigger in 

order to insure that the time constant is reasonable. We hope to describe some practical 
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approches to choosing the arbitrary functions which appear in our family of control laws 
so that the resulting controler will meet specific design criteria in a future paper. 



For a first rough test of our control law, we ran numerical simulations of the system 
using both our control law and a linear control law with the arbitrary functions and 
parameters specified above. We used 10000 diferent sets of initial conditions. It appears 
that the linear controler has a shorter settling time than the nonlinear control law, when 
it does not blow up. However, the nonlinear control law appears to stabilize the system 
starting from any of the initial conditions stabilized by the linear control law, and also 
appears to stabilize the system starting from many initial conditions for which the linear 
control law produces an unbounded response. 




linear controller nonlinear controller 

Figure 2 



Figure 2 shows graphs of the angular position versus time produced using the linear 
control law and the nonlinear control law starting from initial conditions: 6 = .5, 9 = 
— .5, x = 0, and x = 0. The graphs of the other state variables and control input are 
qualitatively similar. This output is typical of the responce obtained when both control 
laws appear to stabilize the system. 
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linear controller nonlinear controller 

Figure 3 

Figure 3 shows graphs of the angular position versus time produced using the linear 
control law and the nonlinear control law starting from initial conditions: 9 = 1.25, 
6 = 1.3, x = 0, and x = 0. For the same initial conditions, Figure 4 shows graphs of the 
cart position versus time produced using the linear control law and the nonlinear control 
law. 
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linear controller nonlinear controller 

Figure 4 

For the linear controller, x(t) increases exponentially. For the nonlinear controller, x(t) 
stabilizes to 0, but the system is underdamped, as we discussed previously. The graphs 
of the other state variables and control input are again qualitatively similar. These 
graphs are representative for initial conditions for which the linear control blows up, and 
the nonlinear control stabilizes. All of our numerical simulations were produced using 
MATLAB 5.1. 



Our infinite dimensional family of control laws contains a one dimensional family of 
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control laws, which was found previously. Bloch, Leonard and Marsden have developed 
a method for constructing control laws for mechanical systems with symmetry, [Bloch, 
Lenard, Marsden, 1997]. Their control law will retain the symmetry, so it cannot produce 
a truly Lyapunov stable equilibrium. It will only produce a Lyapunov stable equilibrium 
in shape space. This means that their method will not work for an inverted pendulum 
cart on a hill, say. If the cart is on level ground, the quantity x + Acos(9)6 will be 
conserved. The quantity x + Acos(9)0 is conserved in the closed-loop dynamics of any of 
the control laws developed by Bloch, Leonard and Marsden. It follows that any solution 
which passes through a point with 9 = and i^O will run off to infinity. 

LINEARLY NON-CONTROLLABLE EXAMPLE 

The final example that we will consider is an abstract system which is open-loop unsta- 
ble with the property that the linearized system cannot be stabilized by any control law. 
Our method will generate a control law with a globally asymptotically stable equilibrium. 

The abstract system is: g = dx 2 + dy 2 , P = dy®J^ and V = -fx 4 + 45x 2 y 2 + 32xy 3 . 
In coordinates, the control problem reads: 

x - 6x 3 + 90xy 2 + 32y 3 = 
y + 90x 2 y + 96xy 2 = u. 

This system is, clearly, open-loop (u = 0) unstable. For example, x = (e -1 — y/St) -1 , 
y = is a solution to the system for any positive e. The linearized version of the system 
is: x = 0, y = u. Clearly, u can have no effect on x. 

To apply our method to this system, let XPX = a-J^ + ^-§^, then the A-Equation (9) 
reads: 

^ a o ^ a 

dx dy 

Thus, a is any constant and \x is any function. Pick a = 1 and \i = 1. The (/-Equation 
(10) becomes: 

dfju , 9gn = Q 
dx dy 

Pick <7n = 2. The equations 'gX = g lead to ^22 = 1 and ^12 = — 1. It is easy to check 
that the model mass matrix is positive definite. The flow equation for V is 

dV dV o , o 

— + — = -6x s + 90xy 2 + 32y 3 . 
dx dy 
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It is not hard to check that 



V = (x 2 - 3xy) 2 + (x 2 - 4xy - 2y 2 ) 2 

is a positive definite solution to this equation. Finally, pick c = (y solution 
to Pic) = P(c) = 0. As before, we show that there is no non-trivial solution to the 
controlled equations with y — x = 0. Thus, (0, 0) is a globally asymptotically stable 
equilibrium. 
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